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UNIT I: FIRST-ORDER DIFFERENTIAL EQUATIONS 



We set forth fundamental principles in the analysis of differential equations. We illustrate the 
use of integration to find the solutions of first-order linear differential equations and special classes 
of first-order nonlinear differential equations, called separable equations. Substitution techniques 
are used in studying linear fractional equations and special kind of second-order differential equa- 
tions. 



LECTURE 1. INTEGRATION AND SOLUTIONS 



We recall the fundamental theorem of calculus 

(1.1) ^iy^'^'^'=^^''^^ 

if / is continuous on an interval xq ^ I. A solution of the differential equation 

is the function y = (f){x) which satisfies the differential equation on /. Upon inspection of (1.1), 
then, y = fis)ds is a solution of (1.2). This leads to an existence result. 

J Xn 



Theorem 1.1. If f{x) is continuous on an internval Xo & I then given an arbitrary number yo there exists 
a unique solution of (1.2) satisfying y{xo) = yo- The solution is given as 

y{x) = yo+ f{s)ds. 

Jxo 

Exercise. Prove the uniqueness. 

Remark. 1. The theorem specifies the interval of existence (xq G /) and the class of functions 
considered (the class of continuous functions). It asserts the existence and uniqueness of a solution, 
prescribed the initial condition y(xo) = yo- 

2. In the statement of the theorem, the interval of existence is /, regardless of the initial condi- 
tion. It is a special property of linear equations. For nonlinear equations, in general, the interval 
of existence depends on the initial value, e.g. the solution of the initial value problem 

^ = y^ y{0) = yo ^ 0, 

is given as y{x) = . It is defined on x e [0, 1/yo) for yo > 0. 

V / yo ) 

rx 

3. The definite integral / f{s)ds is defined as a limit of Riemann sums, as long as / is continu- 

Jxo 

ous; it doesn't need to find a formal expression for the indefinite integral f f{s)ds to give meaning 
to the definite integral, e.g. the error function erf(x) = e~'^^ds and the sine integral function 
Si(x) = /o^(sin s)/sds are commonly defined as definite integrals. 
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As an illustration, the solution of the initial value problem 



= sinx^, y(0) = 

dx 

is given by the Fresnel sine integral function S(x) = sin ds. There is no elementary function F 
such that F'{x) = sin x^, but the function S(x) defined as a definite integral gives a perfectly good 
function. 

The preceding discussion leads to how to solve differential equations of the form (1.2) by in- 
spection. For any xq, one solution is the function / f{s)ds. Other solutions are, then, obtained 

J xo 

by adding an arbitrary constant C to this "particular" solution. Thus, the solutions of y' = 
are the functions y = J ds = {^/7: /2)erf{x) + C. From any one solution curve of (1.2), the 
others are obtained by the vertical translations (x, y) ^ {x,y + C) and they form a one-parameter 
family of curves, one for each value of the parameter C. 



Quadrature. When the solution of a differential equation is expressed by a formula involving one 
or more integrations, it is said that the equation is solvable hy quadrature. The term "quadrature" 
has its historical origin in the connection of integration with area. In plane geometry, a problem 
of quadrature, such as quadrature of the circle is a problem about the area of a plane figure. Not all 
differential equations can be solved by quadrature. In the following lecture, we will show that the 
first-order linear equation 

y +p{x)y = q{x) 

can be solved by quadrature. But, the second-order differential equation 

y" + p{x)y' + q{x)y = r(x) 
cannot be solved, in general, by quadrature, except for some special cases. 



The next simples type of differential equation is 
(1.3, 

Such a differential equation is invariant under horizontal translations (x, y) i— > (x + c, y). Geomet- 
rically, it means that any horizontal line is cut by all solution curves at the same angle (such lines 
are called "isoclines"). Therefore, ii y = 0(x) is a solution of (1.3), then so is y = (p{x + c) for any 
c. The differential equation (1.3) can be solved by writing it as dy/g{y) = dx and integrating. 



Example 1.2. Consider 



dy 2 



(1.4) ^ = _ 1. 

dx 

Since — 1 = {y — l){y + 1), the constant functions y = ±1 are particular solutions of (1.4). 
They are called steady states, stationary solutions or equilibria, in the sense that these solutions are 
independent of x. 

Next, if lyl < 1 then y^ — 1 < 0, and follows dy/dx < 0. That is, the solution curve is decreasing. 
On the other hand, if |y| > 1, then dy/dx = — 1 > 0, and the solution curve is increasing. It 
gives us the qualitative behavior of solutions curves of (1.4). 
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Figure 1.1. Qualitative behavior of solutions of y' = — 1. 



Using the partial fractions and separating variables (we will discuss this technique in detail 
later), (1.4) is written as 

, / 1 1 

2ax = ay 



,y - 1 y + 1 
Then, by integration, we obtain 

1 ± e2(^-'=) r tanh 1 , 

This procedure of separating variables "loses" the particular solutions y = ±1, but it gives all 
other solutions. 

Note that if y = is a solution of (1.4) then so is 1/y = l/(/»(a;). 
Exercise. Discuss y' = y^ — y- 
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